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Abstract 


A cubic spline based Galerkin-like method is developed for the identifi- 
cation of a class of hybrid systems which describe the transverse vibration of 
flexible beams with attached tip bodies. The identification problem is formu- 
lated as a least squares fit to data subject to the system dynamics given by a 
coupled system of ordinary and partial differential equations recast as an ab- 
stract evolution equation (AEE) in an appropriate infinite dimensional Hilbert 
space. Projecting the AEE into spline-based subspaces leads naturally to a 
sequence of approximating finite dimensional identification problems. The 
solutions to these problems are shown to exist, are relatively easily computed, 
and are shown to, in some sense, converge to solutions to the original 
identification problem. Numerical results for a variety of examples are 
discussed. 
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INTRODUCTION 


In this paper we develop an approximation scheme for the identification of 
systems describing the planar transverse vibration of beams with attached tip 
bodies. Standard models from the theory of elasticity for the vibration of 
structures of this type involve hybrid systems of coupled partial and ordinary 
differential equations which describe the motion of the beam and tip bodies 
respectively. The approximation scheme is based upon the formulation of the 
identification problem as a least squares fit to data subject to the dynamical 
equations recast as an abstract evolution equation in an infinite dimensional 
Hilbert space. Using a cubic spline based Galerkin method, a sequence of suc- 
cessively higher (but finite) dimensional state approximations are construct- 
ed. This leads naturally to a sequence of approximating identification prob- 
lems, each of which is shown to have a solution that can readily be computed 
using standard numerical techniques. Results from linear semigroup theory and 
the theory of evolution operators are used to demonstrate convergence of the 
state approximation. This in turn is used to argue that solutions to the fi- 
nite dimensional identification problems, in some sense, approximate solutions 
to the original identification problem. Our effort here is similar in spirit 
to the approach taken in [1], [2], [3], [7], and [11] wherein approximation 
schemes for the estimation of parameters in beam equations with standard boun- 
dary conditions (i.e., clamped, simply supported, free, etc.) are developed. 

Our work is based to a large extent on the ideas suggested in the short note by 
Burns and Cliff [5] . 

Although our general approach is applicable to a broad class of problems 
(see Section 4), to illustrate our method we consider a beam, clamped at one 
end and cantilevered at the other with an attached tip body. In Section 2 the 
derivation of the equations of motion for the beam/tip body system is outlined, 
the equivalent abstract evolution equation is derived and the identification 
problem is formulated. In Section 3 the approximation scheme is constructed 
and convergence results are discussed. Numerical results for several examples 
are presented in Section 4. 

Our notation is, for the most part, standard. The usual Sobolev spaces of 
real-valued functions on the interval [a, b] whose k^ derivatives are 
are denoted by H^(a, b) . These spaces are assumed to be endowed with the 
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usual Sobolev inner products <• , • > k and their induced norms j - |k* For z a 
normed linear space with norm j*J z and f : [0, T] -► Z we say that f e 

2 rT | I 2 

L ([0, T] , Z) if J 0 f(t) z dt < <° . Similarly f will be said to be an element 
of C ([0, T],Z) if the map t + f(t) from [0, T] into Z is Jt time continuously 
differentiable on (0, T). Finally for a function of one or more real varia- 
bles, the symbol D Q f (D*f) will be used to denote the 1 st (k th ) derivative of 
f with respect to the independent variable 0 . If f is a function of a single 
variable only, the subscript may be omitted. On occasion, the short-hand nota— 
tion Dgf(0 o ) or Df(0 Q ) will be used in place of D^fjOg or Df|0 Q to denote the 
derivative of f evaluated at 0g. 

2. THE PARTIAL DIFFERENTIAL EQUATION AND BOUNDARY CONDITIONS, THEIR 

ABSTRACT FORMULATION, AND THE IDENTIFICATION PROBLEM 

We consider (see Fig. 2.1) an inextensible beam of length £, having spa- 
tially dependent linear mass density p and flexural stiffness El. The tip body 
is assumed to be of mass m, have mass center at a distance c from the end of 
the beam directed at an angle 6 measured from the extension of the longitudinal 
axis of the beam and having moment of inertia J about its center of mass . 


*v 



Figure 2.1 . 
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Assuming small deformations (|u(t, x)|<< £), using the standard Euler- 
Bernoulli theory (neglecting rotatory inertia and shear deformations) and ele- 
mentary Newtonian mechanics the equation describing the verticle displacement 
u(t, x) of the beam at position x e[0, £] at time t > 0 

pD 2 u = -D 2 EID 2 u + D TD u + f (2.1) 

t X X X X 

is obtained where T(t, x) is the internal tension resulting from loads directed 
parallel to the beam's longitudinal axis and f(t, x) describes effects due to 
lateral or transverse loading and/or rigid-body rotations (see [6] , [14], 

[17] ). 

If we let S denote the shear force and M the bending moment then using 
the standard moment equilibrium equation for a beam under tension 

S = TD u - D M, 

X X 

the basic bending moment-curvature relationship from the Euler-Bernoulli theory 

M = EID 2 u, 
x 

and the equations for the translational motion of the tip body we obtain the 
first boundary condition at x = £ 

2 2 2 
mD u(t, £ ) + me cos 6 D D u(t f £) = g (t) + D EI(£)D u(t, £) 

r t x 1 x x 

- x(t, £ )D^u (t, £ ) , (2.2) 

where g-j describes the net translational effects on the tip body's center of 
mass which result from externally applied lateral loads and moments (see 

[18] ). The second boundary conditon at x = £, derived from the equations for 
the rotational motion of the tip body is given by 

JD 2 D u(t, Z) = - c cos 6 D EI(£)D 2 u(t, Z) - ElU)D 2 u(t, Z ) 

T- X XX X 

- c sin 6 x(t, Z), + g 2 (t), (2.3) 

where g 2 is defined analogously to g 1 with regard to rotational effects 
(see [18]). 
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The boundary condtions at the clamped end, x = 0, are of course, given by 
u(t, 0) = 0 D x u(t, 0) - 0, (2.4) 

while the temporal boundary conditions (initial conditions) are of the form 


u(0, x) = <|>(x) D t u ( 0 , x) = ip(x). 


(2.5) 


Hie equations (2.1) and (2.2) as they are written above are, in fact, non- 
linear. Indeed, the internal tension x(t, x) is the sum of any externally ap- 
plied loads o(t, x) which are directed parallel to the longitudinal axis and 
the axially directed force me sin 6 D^^uft, £) which results from the angular 
acceleration of the tip body (see [18]). Discarding the nonlinear terms in 

(2.1) and (2.2) as second order effects and choosing w-|(t) = £), 

W 2 (t * = °t D x U(t ' w 3 * tf x) = D x ll * t ' x) ' and w 4 * tf x * = D t u(t ' x) we rewrite 

(2.1) , (2.2), (2.3), and (2.5) in state space form as 


D t w(t, x) 


aD EI(£)w (t, £) + BEI(£)w (t, l) -i.lt, £) w,(t, 6 )d6 

* J j m 0 3 

-pD x EI(£)w 3 (t, £) -yElUJw (t, l) 

°x W 4 (t ' X) v 

-- D x EIw 3 (t, x) +- D x a(t, x) J* w 3 (t, 0 )d0 


pa(t, £) + g (t) - 3g (t) 
ml 2 

Xa(t, £) + Y9 2 (t) 


( 2 . 6 ) 


~ f(t, x) 
P 
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w(0, x) = 


*U> 


D i|>U> 
x 


D 2 <\> (x) 

X 


4^(x) 


where 


a = 


J + me 


2 2 2 ' 
mJ + m c sin 6 


B = 


c cos 6 


2 2 * 
J + me sin 6 


T « 


2 2 9 
J + me sin 6 


M = " 


me cos 6 sin 5 

2 2 f 
J + me sin 6 


X = - 


c sin 6 


2 2 ' 
J + me sin 6 


and 


w(t, x) = (w (t), w 2 (t), w 3 (t, x), w 4 (t, x)f. 


Recalling (2.4) displacement, u(t, x) is recovered from w(t, 


u (t, x) = J* J® w 3 ( t , t) dt d0 . 


(2.7) 


x) by 


The identification problem which we shall consider involves the estimation 
of the flexural stiffness El, the mass density p, the externally applied forces 
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and moments in the form of c, f, g-j , g 2 , and the initial conditions 4> and 
^ • Although (laying identifiability questions aside) our approximation and 
convergence results would be applicable to inverse problems involving the esti- 
mation of any or all of the parameters in (2.6) and (2.7), for ease of exposi- 
tion, we assume that the rigid-body mass properties m, J, c, and 6 of the tip 

body are known a priori. Hie identification problem is formulated as a least- 

squares fit to data. Our approach is based upon recasting (2.6) and (2.7) in 
terms of an abstract evolution equation. 

Let Q be a subset of R^ and assume that the unknown temporally and/or 
spatially varying functions El, p, o, f, g u g 2 , $, and ip appearing in 
(2.6) and (2.7) which are to be identified have been parameterized by qeQ 
(i.e., EI(x) = EI(x; q), p(x) = p(x; q), o(t, x) = a(t, x; q), etc.). We 
require and assume throughout that the following assumptions hold: 

A1 : Q is a compact subset of R L . 

A2 : The mappings q EI(q) and q p(q) are continuous from Q into 

H 2 (0, Z) and Q into H^O, Z), respectively, and there exist posi- 
tive constants m^, m^ , H^, M p such that _< EI(q) <_ 

m < p (q) < M for all qeQ. 

P p 

A3: The mappings q <|> (q ) and q i|;(q) are continuous from Q into 

H 2 (0, Z) and Q into H®(0, Z) , respectively. 

A4: There exists a T > 0 such that the mapping t o(t, •; q) is an 

element of ( [0, T] , H^(0, Z)) and the mapping q <j(t, •; q) is 
continuous from Q into fOfO, Z) for each t£ [0, T] . 

A5: The function f satisfies: 

(i) The mapping t f(t, •; q) is an element of L 2 ( [0, T] , 

H°(0, Z)) for each qeQ. 

(ii) The mapping q f(t, • ; q) is continuous from Q into 
H®(0, Z) for each te [0, T] . 

(iii) There exists Kf£L 2 (0, T) independent of qeQ for which 
| f ( t, •; q)|o <. Kf(t) for all qeQ and te[0, T] . 
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A6s The functions gj_, i » 1, 2 satisfy 

(i) 91^1,2(0, T) for each qeQ. 

(ii) The mappings q gjjt; q) are continuous from Q into R for 
each te[0, T] . 

(iii) There exist Kg^eL 2 (0, T) independent of qeQ for which 
jgi(t; q) j ± Kg i (t) for all qeQ and te[0, T] . 

Let Z = R 2 x H°(0, Z) x H°(0, Z) and for each qeQ let Z q denote 
the Hilbert space {z, <*, *> q } where 

<( V V Vl ), (r 2 r u 2 , v 2 )> q = r^Wr 2 + /J EKqJu^ + J* rfqJv^ 

with 


W = 


m 


me cos <5 


me cos 6 


J + me 


The definition of <•, *>q is motivated by an energy expression. Indeed the 
sum of the kinetic and strain energies for the system described by (2.6) and 
(2.7) is given by 


1 I 1 I 1 /•£ 2 1 ,£ 2 

j [w 1 (t) , w 2 (t)]W I + jl o pw 4 (t ' x) dx + \ Iq EIWj (t f x) dx 

(see [18]). 

Define Ag(q) : D 0 C.Zq ♦ Z q by 

T 

D q = {(r, u, v) ez: u, veH 2 (0, l) , v(0) ■ Dv(0) * 0, r ■ (v(£), Dv(Z)) }, 

A Q (q)((v(A), DvU)) T , u, v) » ( ( aDEI ( Z j q)u(£) + 0EX ( Z; q)u(*), 

- 6DEI (Z; q) u( A) - yEI(*; q)u(£)) T , D 2 v, 


_J 

P(q) 


D 2 EI (q) u) . 


( 2 . 8 ) 
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For each te[0, T] define B(t? q) : Zq Zq by 

B(t; q) (r, u, v) » ( (- ^ o(t, it q) J* u, 0) T , 0, D x c(q) /*u) 

and A(tr q) ; D Q CZ q -► Z g by 

A(t; q) = A Q (q) + B(t; q) . 

For each qeQ it can be argued that the operator Ag(q) is densely defined 
and dissipative (in fact conservative, i.e. < Ag(q) z,z > q = 0, zeDg) . 
Moreover, it can be shown that it is skew self adjoint (i.e., Ag(q)* = 

-Ag(q)) and therefore that it is closed and maximal dissipative. This in 
turn implies that Ag(q) is the infinitesimal generator of a Cg semigroup of 
contractions {Sg(t; q) :t o} on Zq (see [10], Theorems 4.4 and 4.5). It 

is in fact the case that Stone's theorem ([20], pg. 345) implies that Sg(t; 

q) is defined for t < 0 and that {Sg(t; q) < t < is a Cg group of 
unitary operators on Zq. 

The operators B(t; q) are bounded (uniformly in t and q for te(0, T) and 
qeQ) from which it follows that {A(t; q) }te[0,T] is a stable family of infin- 
itesimal generators of Cg semigroups {s^fT; q) :t o} on Zq with stabil- 
ity constants 1 and K = sup |B(t; q) | q (see [12], Section 5.2). Since Dg 

te [0,T] 
qeQ 

is independent of t and t -*■ o(t, •, q)eC^((0, T) , H^(0, £)) for each qeQ, 
the homogeneous initial value problem 

D fc z(t) = A(t ; q)z(t) 0 _< S < t _< T (2.9) 

z(s) - z Q (2.10) 

with zgez has a unique evolution system {u(t, s; q):0<^s<^t<^T} associated 
with it which satisfies 

(i) |u(t, sj q) | q <_ eKtt-s) 


(ii) IJ(t, s; q)D 0 CDg 
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(iii) U(t, s; q)z is strongly continuously differentiable in Z q for all 
zcDq with 

D t U(t, S; q)z = A(t; q)U(t, s; q)z 
and 

D s U(t, S; q)z = -U(t, S; q)A(s; q)z 
for 0 <_ s t _< T. 

If ZqGDq, z(t) = U ( t, s? q)Zg is the unique solution to (2.9) and (2.10). 

For each te[0, T] and qeQ let F(t; q)eZ and z 0 (q)eZ be given by 

F(t; q) = [(po(t, % ; q) + — g (t; q) - Bg o (t; q) , Xo(t, ij q) 

m i z 

+ Y 9 2 (t; q)) T , 0, - p ^ f(t, •; q)) 

and 


z Q (q) = (r (q) , D 2 <{>(q), ip (q) ) 

respectively where r(q)eR 2 and consider 

D^_Z ( t) = A(t;q)z(t) + F(t; q) (2.11) 

z (0) = Z Q (q) . (2.12) 

Writing formally 

z(t; q) = U(t, 0; q)z Q (q) + U(t, T; q)F(x; q) dx , (2.13) 

assumptions A3-A6 imply that the function t z(t; q) is well defined and con- 
tinuous from (0, T] into Zq. If, in addition, t -*■ F(t; qJeC 1 ( [0, T] , Z q ) , 

tj>(q)eH 2 (0, l), 4i(0; q) = Dip ( 0 , q) = 0 and r(q) = (xp (^ ; q) , Dip (X- ? q) ) T 

(i.e., z 0 (q)ED 0 ) then z(t; q) as given by (2.13) is the unique classical 
solution to (2.11), (2.12) in the sense that z(* ; q) eC 1 ( [0, T] , Z q ), 

z(t; q)eD 0 , 0 ^ t _< T and (2.11), (2.12) is satisfied in Z q (see [12], Sec- 

tion 5.5). Under assumptions A3-A6 only, however, a classical solution to 
(2.11), (2.12) does not, in general, exist. In this case, z(t; q) as given by 
(2.13) is known as a mild or generalized solution to (2.11), (2.12) in that it 
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is the limit of classical solutions to sequences of problems of the form 
(2.11), (2.12) for which a unique classical solution does exist. (See [12]). 

In light of the above remarks we use (2.13) to formulate the identifi- 
cation problem. For each xe[0, A] define the operators C(x):Z R by 

C(x) (r , u, v) = J* Jg u(o) da dT. (2.14) 

We assume that we have been provided with displacement measurements , 

{u(t , x )} , t e[0, T], i = 1, 2, ... v, x. e [0, £) j =1,2, ... y, taken 

i J i— I f V 1 ] 

j=i ,y 

from the actual system and state the identification problem as 
(ID) Find qeQ which minimizes 

v y 2 

J(q) = E E jc(x. )z(t. ; q) - u(t., x.)j 
i=1 j=1 1 3 1 1 3 I 

where z(t; q) is given by (2.13). 

Ihe infinite dimensionality of the constraints, (2.13), of course necessi- 
tates the use of some form of approximation in solving problem (ID). We devel- 
op one such scheme in the next section. 

3. APPROXIMATION AND CONVERGENCE RESULTS 

Our approximation scheme is based upon the formulation of a sequence of 
approximating identification problems in which the underlying state equations 
are finite dimensional semi -discrete approximations to (2.13). The approxi- 
mating evolution equations are constructed using a standard cubic spline based 
Galerkin approach to effect the spatial discretization. It will be shown that 
each of the approximating identification problems has a solution, and via con- 
vergence of the states, that the resulting sequence of solutions admits a sub- 
sequence which converges to a solution to problem (ID). 

Working abstractly at first, for each N = 1, 2, ... and each qeQ let Z N 

N q 

be a finite dimensional subspace of Z which is contained in D^. Let P denote 

q 0 q 
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N 

the orthogonal projection of Z ^ onto with respect to the <• , • >^ inner prod- 

N N N N N N 

uct. Define the linear operators A^(q): Z + Z , B (t; q):Z *> Z , and 

0 q q q q 

A N (t? q):Z N + Z N by A*?(q) = P**A (q), B N (t; q) = P N B(t; q), and A N (t; q) = 
q q 0 q 0 q 

N N N N 

P^A(t; q) = A Q (q) + B (t; q) , respectively. Hie finite dimensionality of Z^ 

implies of course, that each of these operators is bounded, although not neces- 
sarily uniformly in N. 

Since A^fq) is conservative, the B(t; q) are bounded uniformly for 

N 

te[0, T] and qsQ and the P^ are orthogonal projections, it follows that the 

N N 

Aq ( q) are conservative and that the B (t; q) are bounded uniformly in N as 

well. Indeed for z N eZ N we have 

q 

N, . N N N , . N N ,/vNN 

<A (q)z , z > = <P A (q)z , z > = <A (q)z , z > = 0 

0 q q 0 q 0 q 

and 


N M M M M 

B (t; q) z = P B(t; q)z < B(t; q)z < K z 

q q q — q — q- 


N 

This in turn implies that the A Q (q) are inf initestimal generators of C Q semi- 

r N 1 N 

groups of contractions, {S Q (t; q):t 2. Oj on Z^ and that the initial value 
problems 


N N N 

D t Z N (t) = A N (t; q)z N (t) 


N,_ N 

z (0) = z Q 


(3.1) 

(3.2) 


r N i 

have unique evolution systems {U (t, s; q) :0 s <_ t <_ T} associated with them 
which satisfy 


(i) 

luV, S; q) 1 

. K(t-s) 
< e 


! iq - 

(ii) 

D t U N (t, s; q)z N 

* t .N a . N . , , N _N 

= A (t; q)U (t, s; q)z for all z eZ , 

q 


0 < s < t < T 



12 


The two parameter families U N (t, s; q) are the solution operators for 
the initial value problems (3.1) and (3.2). We note that since for each N, 

Z q dimensional, once a suitable basis has been chosen, the initial 

value problem (3.1) and (3.2) can be written in matrix form with U N (t, s; q) 

then being represetned by the corresponding principal fundamental matrix solu- 
tion. 

For each qeQ and N = 1, 2, ... we define the function z N (»; q):[0, T] z N 


z N (t; q) 


= U N (t, 0; q)P^z 0 (q) + /Jj U N (t, T; q)p N F( T ; q) dT 


and state the approximating identification problem as 


(3.3) 


(IDN) Find qeQ which minimizes 
N v U 

J (q) = E £ |c(x )z N (t ; q) - u ( t . , x . ) j 2 

i=1 j=1 1 3 1 3 I 

where z N (t; q) is given by (3.3) and the operators C(x) are as they were de 
fined in (2.14). 


Once a basis for Z^ has been chosen, problem (IDN) takes the form of a 
least squares minimization problem subject to a linear non-autonomous, non- 
homogeneous matrix ordinary differential equation which can be solved (assuming 
for the moment that a solution exists) using standard techniques and readily 
available software. 

In terms of the abstract formulation above, our general convergence re- 
sults are summarized in the following two theorems. 

Theorem 3.1 

Suppose {q N } is a sequence in Q with q N q*eQ as N ■*• ». Suppose 
further that 
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N 

(1) P + I strongly in Z uniformly m q for qeQ as S + *. 

q q 

(2) lim |u N (t, s; q N )P N „z - P N „U(t, s; q*)zl = 0 uniformly in t, s 

..I N N N 

N-*-°o 1 q q q 

for 0 _< s <_ t _< T. 

Ihen lim |z N (t; q N ) - z(t; q* ) I -► 0 for each te [0, T] . 

N-v® 1 q N 


Proof 


jz N (t; q N ) - z(t; q* ) I <_ lu N (t, 0; q N )P N N ( z Q (q N ) - z Q (q*))j N 

1 q q q 

+ ju N (t, 0; q N )P N N Z Q (q*) - P^Uft, 0; q*)z Q (q*)j ^ 

1 q q q 

+ |(P N n - I)U(t, 0; q* ) z q ( q* ) I 

q q 

+ Jg ju N (t, t; q N )P N N | N |f(t; q N ) - F(x; q* ) j dT 


q q 


+ /q |(u N (t, X; q N )P N N - pV( fc , T; q* ) ) F ( X ; q* ) j dx 

q q 


+ | 0 j(p N - I) u(t, x; q* ) F ( x ; q* ) I N dx 

q q 


N N 

The properties of U (t f s; q), the fact that is an orthogonal projection 
and assumption A3 imply that the first term above tends to zero as N 00 • 
Hypotheses (1) and (2) in the statement of the theorem imply that the third and 
second terms respectively tend to zero as N -► 00 • Similar arguments and the 
Lebesgue dominated convergence theorem can be used to argue that the last three 
terms tend to zero as well and the theorem is proven. 


Theorem 3.2 

Suppose that hypotheses (1) and (2) of Theorem 3.1 hold. Suppose further 

— N 

that for each N = 1, 2 , ... problem (IDN) has a solution denoted by q . Then 
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H. N, 

the sequence {q } has a convergent subsequence, {q } with q qe Q k ® . 
Moreover, q is a solution to problem (ID) . 


Proof 

_N 

The existence of the convergent subsequence {q } is an immediate con- 

_N 

sequence of Assumption A1 . Theorem 3.1, therefore, implies that |z N (t; q *) - 

z(t; q)| + 0 as k “ for each te [0, T] . lhis in turn implies that 

k 

q 

N N, 

C(x)z (t; q ) - C(x)z(t; q) I -*■ 0 as k «® for each xe [0, Z ] and each 
te [0, T] . It then follows that for any qeQ 


J(q) 


\ N 

lim J (q ) _< lim J K (q) = J(q) 
k-*-°° k-*-°° 


and consequently that q is a solution to problem (ID). The final equality in 
the expression above follows from an application of Theorem 3.1 with the con- 
stant sequence {q}. 

Remark 

In the identification problem (ID) as stated in the previous section, the 
fit is based upon spatially sampled displacement measurements. We note, how- 
ever, that the convergence results given in Theorems 3.1 and 3.2 remain valid 
for identification based upon spatially sampled slope measurements, spatially 
distributed displacement, slope or velocity observations, velocity data at 
x = Z or any combination thereof. 

We next describe a particular realization involving cubic spline functions 
of the abstract ideas presented above and show that the resulting approximation 
system satisfies the hypotheses of Theorems 3.1 and 3.2 

r Ni N+1 

For each N = 1, 2, ... let denote the standard cubic B-splines 

on the interval [0, A] corresponding to the partition A N = {o, — , — , ... Jl} 

N N J 
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(see 

(13)) 

and 

r A Ni N+1 

let { B . ] . „ denote the modified cubic 
1 3 3=1 

B-splines which satisfy 

b"(o> 

_! IN 
= DB. 
3 

(0) 

r "Ni N+1 

= 0, j = 1, 2, ... N+1. The IB.}. 

J i j>2= 1 

are given by 


IN, ^ 
B 1 (x) 

= 

N N N 

bJJ(x) - 2 b"(x) - 2 B^ (x) 



B*(x) 

= 

B^(x) j = 2, 3, ... N+1. 


Let S 

3 <a n ) 

= SPAN S 3 (A N ) = SPAN Z 

N = { (sU ), Ds(£)) T , s, 


a 3 N a a 3 N i N r N i 

s)eZ:seS (A ), seS (A )}, and Z = {z , <• , •> }. Defining 


N r T N h 

t. = ( (0, 0) , B., 0) j = -1, 2, ..., N+1 

D J 


(3.4) 


'T = ((b"(£), Db"(£))\ 0, B\|) j = 1, 2, N+1 

, , N r . Ni N+1 r „Ni N+1 f.NiN+1 ,,r m MiN+1 

we have that Z = SPAN {<!> ^ } j == _ 1 + SPAN l Y jJj =1 ' j J j = _ 1 ' J l^ j J j = i • 

N A q 

is a basis for Z , and since S (A )CS (A )CH (0, £), that Z is a 2N + 4 

q 

dimensional subspace of Z^ which is contained in Dq- 

N . N 

The vector representation £ with respect to the basis (3.4) for P^z, 

where z is an arbitrary element in Zq can be computed using the standard 

N 

normal equation characterization for P : 

q 


N N „ N N 

<P z - z, z > =0 zeZ. 

q q q 


For z = ( r , u, v)e Z a we find that 


. CM N ] - 1 h"(z) 

q q 


(3.5) 
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where 


M 


N 



c N 

q 

0 


0 

IF 


__ 

q 


and 


<<■> 


u 


N 


with 


[c N ]. . 

L q J i+2,j+2 


i; El(q)B^ i , j 


= -1, 0, 1, . .., N+1 


[ D q ] i ^ j = [ ( £ ) , DB^(£)] W 


B N U) 

J 

" N 

DBj (£ ) 


r & N ^ N 

+ J 0 p ( q ) B i B ^ i, j = 1, 2, N+1 




f£ . , N 

J 0 EI(q)uB^ l = -1, 0, 1, ... N+1 


and 


[v N ], = 

1 q J l 


T 

r W 


BAD 

. N 
DB i (£ ) 


+ J o p(q)v^; 


i = 1, 2, ..., N+1 


The matrix representation A^(q) for the operator A^(q) can be computed using 
(3.5) with z = A (q)z N f z N an arbitary element in Z N . We find 

u q 


rtf, x 
A 0 (q) 


[mY 1 k n 

L q J q 
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where 



Similarly the matrix representation B^(t; q) for the operator B N (t? q) is 
found to be 

B^t; q) = [m”]" 1 L*(t) 

where 



with 

[cJJ(t)]. ... = - Jn 0(t ' In B i> 1 = 1 » 2, N+1 

q 1 ' 3+2 ° 1 ° 3 j = -1, 0, 1, N+1 

It then immediately follows that the matrix representtion for A N (t; q) is 
given by 
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If for each te[0, T] and q eQ we set F (t; q) = 

N N -1 N ^ 

w Q (q) = [M q ] H^(z Q (q)), with respect to the cubic 

equation (3.3), in differentiated form, is given by 

sional initial value problem 


Hg(F(t; q) ) and 
splines, the evolution 
the following 2N + 4 dimen- 


D t W N (t) = A**(t; q) W N (t) + F N (t; q) 


(3.6) 


N, . 

w (o) = w Q (q)« 


(3.7) 


The approximating identification problems take the form 
Find qeQ which minimizes 
v y N+1 

i 

j = 

where w^(t? q) is the solution to (3.6), (3.7) corresponding to qeQ and 


J N (q) = 


V \1 N+l 

Z Z I E W V + 2 (t i ; 3)s??(x ) - u(t , x 4>| 

i=1 j = l 'k — 1 k+2 1 * 3 1 3 1 


s^(x) = C(x) b\) do di j = -1, 0, 1, ..., N+1. 

In order to demonstrate that the scheme described above satisfies the 
hypotheses of Theorems 3.1 and 3.2 the following approximation theoretic re- 
sults for cubic splines will be required (see [15], Chapter 4). Let I N de- 
note the standard cubic spline interpolation operator on [0, £] corresponding 
to the partition A N . That is, for <f> a function defined on the interval 
[0, SL ] , I N (j) is defined to be that element in S 3 (A N ) which satisfies 
d N< t> ) (^p) = <l>(-j|p)f j = 0, 1, 2, ..., N and D(I N <f>)(ip) = D<}> ) j = 0 and N. 

For <J>eH 2 (0, l) 


Proposition 3.1 
k , _N 


D R (A - <^)|o 1 C k N ’ 2+k | I)2, < , |o 


k = 0, 1 


.1 • 


where is independent of <j> and N. 
Proposition 3.2 For <|>eH 4 (0, £) 

K' 1 "* - +’|o ± c k N ’ 4+k | D4 ' |, |o 


k = 0, 1, 2 
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2 

where is independent of 4> and N. 
Let 


O n 

V = R x H (0, £) and let V = {v, <<•,•>>} , 

q q' 


« (r i ' V' ‘V V ”, * r i" r 2 + <V 1' Vp(q) 


where 


<*» ^ > p (q) = /o 

Let 

V N = {(n, <)>)eV:<|>eS 3 (A N ), n = (<{.(£), Dcf> { A ) ) T } 

and 


N r N l 

v q « (v , «•, ->> q }. 


N 0 3 N 

Let (q) denote the orthogonal projection of H (0, Z) onto S (A ) with respect 

rZ N 

to the inner product <<j>, = Jq and P 2 ^q) denote the orothogo- 

nal projection of V on to V N . Adopting the convention that for z = 

N N N N N 

(r, u, v)eZ , z = P z will be denoted by (r , u , v ), it is easily seen 

q q 

. N N , , , , N N . N , w x 

that u = P^ (q)u and (r , v ) = P^fqHr, v) • 


Lemma 1 Let z = (r, u, v)eD 0 . Then 

(1) |u N - uj 2 + 0 as N + 00 uniformly in q for qsQ. 

(2) J v N - vj 2 *► 0 as N + 00 uniformly in q for qeQ. 

Proof 

To verify (1) we show that Jd^(u N - u) j ^ + 0 as N 00 , k = 0, 1, 2 
uniformly in q. Recalling assumption A2, the fact that zsDq implies 
ueH 2 (0, Z) and Proposition 3.1 above we have 
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K • u lo ■ K lq,u - u lo s. K< q )“ - “| EI(q , 


I 1 ” 11 • °U> ^ I 1 "" - “lo 


El 1 -2 | 2 | 

< C N Id u 0 as N + », 

“ m EI 0 1 10 

Hie convergence of the first derivative is argued using the Schmidt inequality 
([15] Theorem 1.5), 


D(u N - u) £ d(u N - I N u + D(I N u - u) 

v u 0 


“* N N 1 -1 o 

£ KNu - I u Q + C^N D 2 u 


£ Kn|u N - u| Q + KN|l N u - uj Q + C^N -1 |d 2 u| q 


£ ^ c 0 1 »* 1 |'> 2 u| 0 + K c’n- 1 |d 2 u| 0 + c;h-’|d 2 u| 0 . 0 as H . 

b I 

For the second derivative, we first note that for ueH 2 (0, £), the Schmidt 
inequality and the first integral relation ([15] pg. 52) imply 

) (q )u | q £ 2 |d 2 (P^(q)u - I N u)| 2 + 2 |d 2 I N u| 2 

£ 2KN |p”(q)u - I N u| 2 + 2 |d 2 u| 2 - 2|d 2 (u - I N u) | 2 
£ 2KN 4 |pN(q) u _ u | 2 + 2 KN 4 |l N u - u| 2 + 2|d 2 u| 2 


£ k 2 D 2 U 2 


(3.8) 


where K is independent of u, N, and qeQ. Now for weH 4 (0, £) 


with w N = P^(q)w, 
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Proposition 3.2 together with the Schmidt inequality imply that |d^(w N - w) = 
— 4+k 

0(N ) , k = 0, 1, 2. Therefore 

|d 2 (u N - u)j 0 £ |d 2 (u N - w N )j 0 + |d 2 (w N - w)j Q + |d 2 (w - u)| Q 

<_ |d 2 (w N - w) j + (1 + k) |d 2 (w - u) j 0 

£ 0(N -2 ) + (1 + K) | D 2 (w - u)| 0 

1 2 N N | 4 

D (u - w ) . Since H (0, A ) is dense in 

H (0, £ ) we can choose w and then N (since the 0(N ) term depends upon D w ) 
to make the right hand side of the last inequality above arbitrarily small. 

Turning next to statement (2) and recalling that zeD q implies that 
veH 2 (0, l), v( 0) = Dv ( 0 ) = 0 and r = (vU), DvU)) T we have 


CM 

JZ 

1 

IV - V . 

< 

1 0 

— m 

r\ 


— m~ ( | C CX N v(£. ) , DI N vU)) T , I N v) - (v(A), Dv(A)) T , 


1 | IN 1 2 1 | _N 1 2 

n>p I v ” v |p(q) m p I v v |p(q) 


< ^ |*" — IS i 

P P 



2 

q 


as N ■> 00 where for <J> a function defined on [o,Jt], I <{) denotes that element in 

s 3 (A N ) which satisfies I $ (^1 ) = <f> (^), j = 1, 2. .N and DI N <|> (A) = D$ (A). 

N N 


The convergence of the derivatives is verified in essentially the same 
manner as it was in the proof of Statement 1 . 

Theorem 3.3 

For the cubic spline based scheme described above, hypothesis (1) of 
Theorem 3.1 is satisfied. 
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Proof 

Let z = (r, u, v) eDq and let z N » P^z = (r N , u N f v N ) . Then r = (v(£) , 
Dv(£)) T , r N ■ (v N (£), Dv N (£)) T , and 


„ i N ,2 

z = | z - z| 

q I q 


= (r N . r) T w(r N _ r) + ju H _ u j 2 


EI(q) 


+ lv N - vl 2 


- " V(£) | 2 + 

u i N |2 i N 1 2 

+ \il u - “lo * M p| v - v lo 


< K,|u N - U |2 ♦ S 2 |vN - ,|| 


where j | w | 1 2 the Euclidean matrix norm of w and K 1 and K are constants 

which are independent of N, z, and qeQ. Lemma 3.1 implies that the right 

hand side of the final inequality above tends toward zero as N -*■ <*> and 
N 

consequently P z z as N ~ for all zeD_. However, D. dense in z (uniform- 
4 0 0 q 

N 

ly with respect to qeQ) and the uniformly bounded in N (being orthogonal 
. . N 

projections) imply that Pz + zasN + » for all zeZ uniformly in q for 

q q 

qeQ and the theorem is proven. 

In order to verify hypothesis (2) of Theorem 3.1 we require the following 
lemma. 

Lemma 3.2 

Let {q N }CQ with q N + q*eQ as N -► Let be the cubic spline based 

N M 

subspaces of Z defined above, P the orthogonal projections of Z onto Z , 

q q q q 

A> = PgA Q (q), and {s^(t; 3 ) 8t 2. o} the C Q semigroups of contractions gen- 
erated by the A^(q) . Then 
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| S 0 (t , qN)P V ■ pN N S 0 (t; qMz | N * 0 

q q q 

as N 00 for each zeZ, uniformly in t for te[0, T] . 

Proof 

Using a variation of the well known Trotter approximation theorem (see 
[4] f Theorem 6.2) the desired conclusion will follow if we can show that 

|r(X, A o (qN))pN N Z " pN N R(X ' A o (q * ))z | N * ° 

q q q 

as N + 00 for each zeZ for some X > 0 where R(X, A) = (X - A) -1 . However 
|s N (t; q) j < 1 implies (see [9]) 


|R(X, A R (q N ))P N N z - P N n R(X, A 0 (q*))z| R 


= R(X # A^(q N ))(A R (q N )P N N - P> 0 (q* ) )R(X , V q *>>*| H 


q q 


< X | (A 0 (q )p N - P N V q * ))R(X ' A 0 (q * ))z | H 
q q q 


N N N N 

" | (A 0 (q )P N- P N A 0 (q * ))y | N 


q q 


where we have chosen X = 1 and y = R(X, A o (q*)zeD Q . Now 


I /* N , N. N N , I 

,( V q )P N - P »V q ,,y N 

rr /t 1 rw 


q q 


. N , N, N N . , 

* <f (A 0 ‘ q > P » - P N A 0 (q lly | N 


q q 


< |(A 0 (q N ) - A 0 (q*))p"„y| „ + |a 0 (q*)(P N „ - l)y| „ 


q q 


N N 

■ T 1 + V 


Recalling the definition of the operator Ag(q), (2.8), the fact that ye Do 

N 

and the estimates given in Lemma 3.1, standard estimates yield 0 as N -*• 09 
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N N 

while q + q* as N + 00 and assumption A2 imply + 0 as N + “• 

Theorem 3.4 

Hypothesis (2) of Theorem 3.1 is satisfied by the cubic spline scheme. 

Proof 

Since for zez and 0_<s_<t_<Twe have 

D^_U( t f s; q)z = A(t? q)U(t, s; q)z 

= Aq (q ) U( t / s; q) z + B(t; q)U(t, s? q)z 


it follows that 


U(t f s? q) z = S Q (t - S; q)z + J* S Q (t - T; q)B(x; q)U(x, S; q)z dr. 
Similarly 

N NN N/-tN MM M 

U (t, s? q)P z = S (t - s? q)P z + / S (t - X; q)B (x;q)U (x , S; q)P z dx . 

q 0 q s 0 q 


Therefore, letting 


A N (t, s) = 


__N m . N • N N * , » 

U (t, S; q )P Z - P U(t, s; q*)z , 
N N N 


we have 


A N (t, s) <_ Isjj (t - s; q N )P N N Z - P^S (t - S; q*)zl 

q q q 


+ /s | s o (t ■ T; qN) | n| bN(t; qN) | n aN(t ' s) dT 
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or 


+ /s | s 0 (t • T? qN) | n| bN(t? qN) | N | (pN N " I)0(T ' S; qMz | N dT 

q q q q 

+ /: | S 0 (t ■ T; qN) | N I pN N I n| (B( t; qN) " B(T? S*))U(x, s ' ( 5* )z | N dT 

q q q q 

+ Jg |(Sq (t - t? q N )P N N - pN N s 0 ( t - X; q* ) )B (t ; q*))U(x, s; q* )z I dx 

q q 'q 

+ Jg | (P N n - DS 0 (t - X; q* )B(x ; q*))U(x, s? q*)z| N dx 

q q 

N l N N N N 1 

A (t, S) £ |s 0 (t - s; q )P N Z - P N S 0 (t - S; q* )z | N 


q q 

+ K f t A N (t, s) dT + / t h N (x; t, s) dx 
' s ' s 


(3.9) 


where K = sup |B(t; q)|^. Lemma 3*2 implies that the first term on the right 
hand side of (3.9) tends toward zero as N + 00 uniformly in s and t for 0 <_ s <_ 
t _< T, while Lemmas 3.1 and 3.2, q N ■* q* as N + ®, assumption A2 and the 
boundedness and strong continuity of the operators imply h N (T; t, s) + 0 as 
N 00 uniforming in t, t, s for 0_<s_<T_<t_<T. Therefore (3.9) can be 
written as 

A (t, s) < e + K I A (t, s) dT 
— 4 s 

with + 0 as N -*• 00 uniformly in t and s. An application of the Gronwall 
inequality yields 


A N (t, 


< e N e K(t - s) 


from which the theorem immediately follows. 

Finally, for the cubic spline scheme, under our general assumptions, using 
either standard continuous dependence results for ordinary differential equa- 
tions or the Trotter approximation theorem it can be argued that for each N = 

1, 2, ... and te [0, T] , z N (t; q) given by (3.3) is continuous in q. Conse- 
quently J N (q) is continuous in q, which together with assumption A1 implies 
that problem (IDN) has a solution. 
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Having now demonstrated that the hypotheses of Theorems 3.1 and 3.2 are 
satisfied by the cubic spline scheme , we turn next to a discussion of examples 
and numerical results which provides an indication of how well the scheme per- 
forms in practice. 

4. NUMERICAL RESULTS AND CONCLUDING REMARKS 

In this section we present numerical results that were obtained by apply- 
ing variations of the cubic spline based scheme discussed in the previous sec- 
tion to the identification of a variety of hybrid systems involving the vibra- 
tion of beams with attached tip bodies. Since our primary objective was to 
demonstrate the feasibility of our scheme, we considered only relatively simple 
examples . 

In the first example we consider a cantilevered beam with a tip (point) 
mass. The second example involves a free-free beam with an attached tip body 
at each end. A cantilevered beam with a tip body subject to an axially direct- 
ed base acceleration is considered in the final example. Strictly speaking, 
the theory developed in the previous two sections applies directly only to the 
last example. However, the relatively minor modifications which are required 
to make our general method and the corresponding convergence results applicable 
to the other two examples should be immediately clear. 

In general the observational data upon which the fits were based was ob- 
tained by generating solutions using fixed (or so called "true") values of the 
parameters with a Galerkin method and a finite number of the unforced system's 
natural mode shapes. Computing the modal frequencies and corresponding mode 
shapes for systems of the type considered here, in general requires the locat- 
ing of zeros of transcendental equations involving the various beam and tip 
body parameters which appear in the problem. The resulting modal equations 
tend to be stiff and must be integrated using an appropriate method if a valid 
solution is to be obtained. 

The approximating identification problems ( IDN) are solved using the IMSL 
(see [8]) routine ZXSSQ. This routine is an interative Levenberg-Marquardt 
Newton's Method/Steepest descent hybrid algorithm for the minimization of the 
sum of squares of a system of functions of several variables. The required 
gradients and entries in the Jacobian matrix are computed numerically using 
finite difference approximations. The method required that we supply initial 
start-up values for the unknown parameters and a subroutine which evaluates 
J N (q) for a given value of q. The latter requirement necessitates the inte- 
gration of the initial value problem (3.6) and (3.7). This is accomplished 
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using the IMSL routine DGEAR , a variable order ADAMS predictor-corrector meth- 
od. It is interesting to note that uinlike the modal approximations, the 
spline equations did not require the use of the stiff option. 'This of course 
makes the spline schemes attractive from a computational point of view. We 
note also that due to the narrow support of the B-splines, the matrices which 
appear in the resulting Galerkin equations tend to be banded, thus facilitating 
efficient integration of the system of differential equations. The inner prod- 
ucts which determine the entries in the mass and stiffness matrices, as well as 
the generalized Fourier coefficients required to project the non-homogeneous 
term and the initial conditions were computed using a composite two point 
Gauss-Legendre quadrature rule. 

Example 4.1 

We consider a cantileverd beam of lenth £, constant stiffness El and mass 
density p with a tip mass of mass m (see Figure 4.1). The differential 
equations and boundary conditions which describe the system are given by 


PD t u = 

-EID U + f 
X 

2 3 

mD. u(t, Z) = EID u(t, Z) + a 

z 

X 

u(t, 0) 

2 

= D u(t, 0) = D u(t, £) 


X X 

u(0, x) 

= 4>(x) D^u ( 0, x) = I)j(x) 


where we have assumed that no axially directed loading is present. 



Figure 4.1 
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We set i * 1,0, El = 1.0, p = 3.0, m = 1.5, assumed that the system was 
initially at rest (i.e., <J> = = 0) and excited the structure with an implusive 

lateral force at the end of the beam at time t = 0. The input disturbance was 
modeled as f(t, x) = 20e“2t e -20(1 -x) an< j g(t) = o. Using the first three 
normal modes of the unforced system to generate observations in the form of 
displacement measurements at positions Xj = 0.125(j+3), j = 1, 2, ..., 5, and 
times tj = 0.2 j, j =1, 2, ..., 10 we identified El and p. The start-up 
values for the Levenberg-Marquardt routine were taken to be EI 0 =0.7 

and p Q = 2.7. The final converged values for El , p and the residual sum of 
squares are given in Table 4.1 below. 


Table 4.1 . 


N 


-N 

P 

0* 

2 

0.9976 

3.0262 

0.39 x 10“ 5 

3 

0.9994 

3.0382 

0.46 x 10 -5 

4 

0.9951 

3.0544 

0.48 x 10 -5 

5 

0.9961 

3.0409 

0.40 x io -5 

6 

0.9995 

2.9976 

0.35 x 10 -5 


True Values: 1.0 3.0 

It is clear that relatively accurate estimates of the parameters can be 
obtained using small values of N. 

Based upon the scheme's performance on examples for which exact solutions 
were available, we feel that the somewhat eratic convergence exhibited in Table 
4.1 is most likely a consequence of using approximate solutions to generate 
observations. 
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Example 4.2 

In this example we identify the stiffness and mass density of a free-free 
beam of length £ which has a tip body (having different mass properties) at- 
tached to each end (see Fig. 4.2). The system was assumed to be initially at 
rest and then excited by a time varying, spatially distributed transverse load 
given by f(t, x) = 10 sin (2TTt)e x . 


y 



Figure 4.2. 


The vibrations of the beam are described by the following partial 
differential equation and boundary conditions (see [16]) 

pD 2 u = -EID 4 u + f t >0 xe(0, £) 

t x 


2 3 2 

Du = -a EID u + 3 EID u 
t 1 x 1 x 

D 2 D U = -B EID 3 u + Y EID 2 U 
t x 1 x 1 x 


t > 0 x = 0 

t > 0 x - 0 


(4.1) 


* ■ 


D 2 D u 
t x 


a_EID 3 u + 3„EID 2 u 
2 x 2 x 

3 2 

= -3„eid U - Y^EID u 
2 x '2 x 


t > 0 x = £ 

t > 0 x = £ 
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u = 0 


V = 0 


t = 0 xe [0, A ] 


with 


a. 


J. + m. c. 
x IX 


2 2 

m, ( J. + m. c. sin 6 . ) 

XX 11 i' 


c. cos 6. 
x x 


2 2 

J. + m.c . sin 6 . 
xxx 1 


and 


Y i 


1 


2 2 

J. + m.c . sin 6 . 
xxx 1 


i = 1, 2 


where the mass properties for the tip bodies, itl, c^ and 6^^ are as they 

were defined for the single tip body problem in Section 2. 

Settling El = 1 .0, p = 3.0, 

m i = 0.75, J i = 0.6, c i = 0.1, 6^ = n/6, 

m 2 = 1,5 » J 2 = °' 4, C 2 = °* 2 ' 6 2 = 1T/3r 

and l = 1 .0 we based our fit on velocity data at the ends of the beam only gen- 
erated at times tj = 0.2 j, j = 1 , 2, ..., 10 using the first six natural 
modes including the two which correspond to rigid body translation and rotation 
(see the remark following the proof of Theorem 3.2). We note that the system 

in (4.1) is free with respect to the intertial frame in which displacement is 

measured. A formula analogous to (2.14) therefore can not be used to recover 
displacement from the solutions to either the abstract or the approximating 
evolution equations. 

The start-up values used were EI 0 =0.7 and p 0 = 2.7. Our results are 
summarized in Table 4.2 
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True Value: 


Table 4.2 


N 

— N 
El 

-N 

P 

a* 

2 

0.9957 

3.0052 

0.26 x 10 -3 

3 

0.9963 

3.0065 

0.18 x 10- 3 

4 

0.9982 

3.0037 

0.10 x 10 -3 

5 

1.0018 

3.0008 

0.71 x 10" 4 


1.0 

3.0 



Example 4.3 

For our final example, we consider a problem of the form discussed in 
Sections 2 and 3. The parameters to be estimated are the spatially invariant 
stiffness El and mass density p. We assume that the system is initially at 
rest and is then acted upon by an impulsive lateral force at the end of the 
beam at time t = 0 and a piecewise constant base acceleration described by 


f(t, x) 

and 


a o (t) - 


20 .-V 01 '-' 

1.0 0 <_ t <_ 1.5 

0 1.5 < t < 3.0 

1.0 3.0 _< t < 4.0 

0 4.0 < t 


respective ly . 

The externally applied axial load a(t, x) can be related to the base 
acceleration ag(t) as follows. Recalling that t denotes the internal tension 
in the beam, equilibrium in the x direction yields D x x = pag, or 


t ( t, x) = -p ( £ - x) a Q (t) + x (t, £). 


(4.2) 


For the tip body we have 


= ma.(t) - me sin 6 D D u(t, 
0 t x 


1 ). 


-T(t, £) 


(4.3) 
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Combining (4.2) and (4.3) we obtain 

T(t, x) = -p(£ - x)a Q (t) - ma Q (t) + me sin 6 • D^D^uft, £), 

from which we immediately conclude that a(t, x) = -a 0 (t)(p(£ - x) + m). 

Taking El = 1.0, p = 2.0, m = 4.0, J = 0.4, c = 0.2, 6 = tt/ 3 and % = 1.0, 
displacement data at positions Xj = 0.75, 0.87, 1.0, at times tj = 4.0, 

4.5, 5.0 was generated using the first 3 natural modes for the unforced, fixed 
base beam/tip body system. The start-up values were chosen as Elg =0.7 and 
Po = 2.5. The results are given in Table 4.3. 




Table 

4.3. 


N 

— N 

— N 

N 

El 

P 

J 

2 

1*0003 

2.0635 

0.35 x 10“ 3 

3 

1.0014 

2.0668 

0.22 x 10 -3 

4 

1.0007 

2.0404 

0.26 x 10~ 3 

5 

.9994 

2.0274 

0.35 x IQ -3 


True Values: 1.0 2.0 

Once again the eractic convergence is most probably attributable to using 
an approximate solution to generate observations. 

We note that strictly speaking, our theory requires that ageC^ . How- 
ever, as is evidenced by the example above the scheme appears to perform satis- 
factorily on problems involving ag which are discontinuous. 

Although it was not considered in the present effort, it is possible to 
include terms which model viscous damping in the partial differential equations 
and boundary conditions. Furthermore, it is possible to extend our approxima- 
tion theory and corresponding convergence results to include the ability to 
estimate parameters associated with the damping effects (see [1], [7], and 
[ 11 ]). 
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It is also possible to develop an approximation theory which is similar in 
spirit to the one presented here for the identification of models based upon 
the use of the higher-order Timoshenko theory to model the vibration of the 
beam. The Timoshenko theory includes effects due to shear deformation and 
rotatory inertia (see [2], [6], [7], and [19]). 

Finally we note that while the scheme described in Section 3 is cubic 
spline based, because of the choice of state variables, it in fact relies upon 
quintics to represent displacement. We are currently developing a method using 
a somewhat more direct approach than the abstract operator formulation employed 
here which does use cubic splines to represent displacement. Starting with a 
weak/variational form of the underlying hybrid system, a cubic spline based 
Galerkin approach is used to construct a sequence of approximating identifica- 
tion problems wherein the constraints are given by a finite dimensional linear 
second order ordinary differential equation in the approximate displacement. 
Based upon preliminary results this scheme promises to be computationally more 
efficient than and to perform as well as the one which was presented above. 
Moreover, it will permit the identification of spatially varying El and p under 
for less stringent smoothness and continuity hypotheses than the ones given in 
assumption A2. This work will be discussed in a forthcoming paper. 
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